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Abstract 

The cosmic microwave background (CMB) is a remarkably distor- 
tionless blackbody, and this strongly constrains the amount of energy 
that can have been injected at high redshift, thereby limiting the role 
that hydrodynamical amplification can have played in cosmic structure 
formation. The current data on primary anisotropies (those calculated 
using linear response theory) provide very strong support for the gravi- 
tational instability theory and encouraging support that the initial fluc- 
tuation spectrum was not far off the scale invariant form that inflation 
(and defect) models prefer. By itself, the (low resolution) 4-year DMR 
data allow relatively precise as normalization factors for density fluc- 
tuation spectra and rough information on the large scale slope of the 
anisotropy power, thereby focusing our attention on a relatively nar- 
row set of viable models. Useful formulae relating the DMR band- 
power to (Tg and post-inflation scalar and tensor power spectra mea- 
sures are given. Smaller angle data {e.g., SP94, SK94) are consistent 
with these models, and will soon be powerful enough to strongly select 
among the possibilities, although there remains much room for surprises. 
In spite of foregrounds, future high resolution experiments should be 
able to allow precise determination of many combinations of the cosmo- 
logical parameters that define large scale structure formation theories: 
mode (adiabatic/isocurvature, gravity wave content), shape functions 
{ns{k),nis{k),nt{k)), amplitudes {a8,C^^ /C2^'^), and various mean en- 
ergy densities {fi, h, ils, f^A, Zreh, fimi/}- Secondary anisotropies arising 
from nonlinear structures will be invaluable probes of shorter-distance 
aspects of structure formation theories. 



in The Evolution of the Universe, pp. 199-223, 
ed. S. Gottlober & G. Borner, (Chichester: Wiley) (1997) 



1 



1 Basics of CMB Anisotropy 



We are in the golden age for cosmic background radiation research, with signals 
unveiled by very high precision spectrum and angular anisotropy experiments 
revealing much about how structure arose in the Hubble patch in which we 
live. The main goal of theoretical anisotropy research is to work out detailed 
predictions within a given cosmic structure formation model of primary and 
secondary CMB temperature fluctuations as a function of scale {e.g., 
Primary anisotropics are those that we can calculate either fully with linear 
perturbation theory, or, as in the case of cosmological defect models, with lin- 
ear response theory of nonlinear seed fluctuations. Because of the linearity, 
primary anisotropies are the simplest to predict and offer the least ambiguous 
glimpse of the underlying fluctuations that define the structure formation the- 
ory. With detailed high precision observations, we expect to be able to use 
CMB anisotropies to measure various cosmological parameters to remarkable 
accuracy (e.g., |, |, |, |, 0, |, |, 0). 

Accompanying spectral distortions to the CMB that may be generated 
during the evolution of nonlinear objects, there will be inevitable secondary 
anisotropies that carry invaluable information about the epochs that the rele- 
vant structures formed. Even if the angle-averaged distortions are well below 
the level that absolute spectrum experiments like COBE's Far Infrared Abso- 
lute Spectrophotometer (FIRAS) probe it is certain that these secondary 
anisotropies are accessible to experiments {e.g., |jl|): the question is only for 
what fraction of the sky do they rise above experimental noise and the primary 
signal. 

To relate observations of anisotropy to theory, statistical measures quite 
familiar from their application to the galaxy distribution have been widely used. 
Denote the radiation pattern as measured here and now by the two-dimensional 
random field AT{q), where —q = {9, (j)) is the unit direction vector on the sky 
(and q is the direction the photons are travelling in). For CMB anisotropies, 
it is natural to expand the radiation pattern in spherical harmonics Yirn{0, 4>) 
and define an ensemble-averaged angular power spectrum: 

AT 

— {q) = ' = l{i + l)(a,*„a,„,)/(27r) . (1) 

At high Ci corresponds to the power in a logarithmic waveband (iln(^). If 
the temperature pattern is statistically isotropic, then {a'^^ai'm') — unless 
i = £' , m = m' . If the initial fluctuations are Gaussian so is the primary CMB, 
hence Ci is all that would be needed to characterize the anisotropy statistics. 
Equivalently, the Gaussian patterns are completely specified by the associated 
2-point correlation function, C{0) = {AT{q)AT{q'))/T^, where cos(6') = q ■ 
q' . (If the statistics are not Gaussian, then an infinite number of A'^-point 
correlation functions are required to specify the statistical distribution.) C{0) 
and more generally the rms temperature anisotropies associated with an £- 
space filter Wg can be expressed in terms of a "logarithmic integral" I[fe] of a 
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Even if an experiment has perfect resolution and all-sky coverage, because 
the observed sky is just one realization from the ensemble the derived Ci and 
C{9) would differ from the ensemble-averaged ones. This effect is called 'cosmic 
variance' and for example implies that, multipole by multipole, the uncertainty 
is AC, = C,/(^+i)i/2. 

Data from an anisotropy experiment are usually expressed in terms of mea- 
surements (Ar/T)p of the anisotropy in the p*'* pixel and a pixel-pixel corre- 
lation matrix Copp' giving the variance about the mean for the measurements. 
The signal {AT/T)p can be expressed in terms of linear filters J-pjm acting 
on the multipole components, aim- {^T/T)^ = J^im^p^rnaem, where Tp,£m 
encodes the experimental beam and the switching or modulation strategy that 
defines the temperature difference. The former filters high £, the latter low £. 
A given theory with power spectrum Cxi has a pixel-pixel correlation matrix 
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We define the band-power of the experiment to be the anisotropy power across 
the average filter Wf. 
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Usually the band-power is the quantity that can be most accurately determined 
from the experimental data. Estimates of band-powers derived for recent ex- 
periments (up to March 1996) are shown in Fig. |l| 

To determine band-powers for an experiment, a local model of Ci is con- 
structed, assumed to be valid over the scale of the experiment's average filter 
Wi. A popular 2-parameter phenomenology has a broad-band tilt j/at as well 
as a broad-band power: 
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As the data improves, a parameterized sequence of best- fit C^'s will be prefer- 
able. 



Because there are so many detections now, Fig. |l| is split into an upper 
and lower panel for clarity, with the upper giving an overview, for experiments 
ranging from dmr at the smallest £ to ovro at the highest and the lower 
panel focussing on the crucial region of the first few peaks in Ci. Data points 
denote the maximum likelihood values for the band-power, the error bars give 
the 16% and 84% Bayesian probability values (corresponding to ±lcr if the 
probability distributions were Gaussian), and upper and lower triangles denote 
95% confidence limits unless otherwise stated. The horizontal location is at 
{tj^ and the horizontal error bars (where present) denote where the filters 
have fallen to e^°'^ of the maximum. The filters W i for the experiments are 
shown in the middle panel. 

For the CMB data sets that have been obtained to date, including COBE, it 
is possible to do complete Bayesian statistical analyses. To determine the best 
error bars on the parameters of a target theory with correlation matrix Cxpp' , 
a recommended method for this analysis llj, 0, |l^ is to expand in signal-to- 
noise eigenmodes, those linear combinations of pixels which diagonalize the 

1/2 1/2 

matrix Cn CxCn , where the noise correlation matrix C„ — Cd + C'res 
consists of the pixel errors Cb and the correlation of any unwanted residuals 
Cres , whether of known origin such as Galactic or extragalactic foregrounds or 
unknown extra residuals within the data. This facilitates the many inversions 
of Cn + Ct required to evaluate the likelihood function, and can also be a 
powerful probe of unknown residuals contaminating the data. The S/N mode 
expansion was used to get most of the bandpowers and their error bars shown 
in Fig. |l|. 

With uniform weighting and all-sky coverage, the S'/iV-modes are just the 
independent Re{aim) and Im^aem)- The uniform noise assumption has been 
used recently to address the ultimate accuracy that satellite experiments might 
achieve [g |l|, |, § |, 0, and is used here in Figs. and | for that purpose. 
The target power spectrum has Cg determined within a la deviation ACe given 

by 

ACt£ ~ I , (D) 

^(^ + \)fsky ^cosh(Aln£)[l + (^-h i)sinh(Aln^)] 

I^Di = — (^u ■ 

Ztt 

If only a fraction fsky of the sky is covered, then for high £, so that the angular 
scale is small compared with the patch probed, the effective pixel number 
scales by fsky ■ The errors are those appropriate for logarithmic binning of width 
±iAln^ about In^, with ^Alnf = 0.05. The filter function associated with 

^This generally introduces smoothing functions but these are nearly unity if A Inf 1. If 
A In £ is so small as to encompass only one £ we recover the usual {£+ ■^)^^^^ cosmic variance 
result. 
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Figure 1: (previous page) Band-power estimates derived for the anisotropy 

data up to March 1996. The lower panel is a closeup of the first two 'Doppler 

peaks'. The theoretical primary power spectra are normalized to the 4-year 

rfmr(53-f90+31)(a-Hb) data. A "standard" ns = l Q = 1 CDM model with 

normal recombination, h = 0.5, fi^h^ = 0.0125 |l9), is the upper solid curve. 

It has Ug = 1.20±0.08. An (almost indistinguishable) dotted curve has the same 

parameters except that it includes a light neutrino with Q„ii, = 0.2 (and i^cdm = 

0.75). It has (78 = 0.83 ±0.06. The upper dashed curve is a vacuum-dominated 

model with a 13Gyr age and h = 0.75 (17a = 0.73, fls ^ 0.02, fl^dm = 0.24). It 

has as = 1.03±0.07. The model whose peak is shifted to high i is an open CDM 

cosmology |0 with the same 13 Gyr age, but now Hq — 60, and iltot = 0.33 

(and Qcdm = 0.30, Qb ^ 0.035). It has as = 0.44 ± 0.04. By Ho = 70, Qtot 

is down to 0.055 at this age. The lower solid curve is a CDM model with 

reionization at z^eh — 30, and almost degenerate with it is a ns— 0.95 standard 

CDM model, whose gravity wave induced component is also shown. Power 

spectra of SZ maps constructed using the peak-patch method |23j are shown 

for a as = 1 standard CDM model and a more realistic (78=0.7 hot/cold hybrid 

( sz\ 

model (solid, Q.mv—Q-i) and a Ua—OS tilted CDM model. Although Cg may 
be small, because the power for such non-Gaussian sources is concentrated in 
hot or cold spots the signal is detectable. Spectra for a dusty primeval galaxy 
model that satisfies the FIRAS constraint is also shown, the larger (arbitrarily 
normalized) part a shot-noise effect for galaxies with dust distributed over 
lOkpc, the smaller a continuous clustering contribution, including a nonlinear 
correction. Average filter functions for a variety of experiments are shown 
in the middle panel. Error curves {la) on the Zreh — 30 model assume a 
homogeneously weighted all-sky survey with a 10' beam and a Q^xK per 10' pixel 
noise, which dominates at high £, while cosmic variance {t) ^ dominates at 
low Ah\i — 0.1 was used. Current band-powers broadly follow inflation- 
based expectations, but could include residual signals from systematic effects 
such as sidelobe contamination. Galactic effects such as bremsstrahlung and 
dust, or secondary anisotropy signals. 
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for a Gaussian beam. It has been divided out to show that the efTective noise 
level picks up enormously above Ig. The parameter Oy is the error-per-pixel 
times dyixei- 

The lowest primary anisotropy curve in Fig. 0, the z^eh = 30 model, has a 
set of one-sigma error bar curves (dotted) on it associated with this uniform 
all-sky coverage error, at small angles due to cosmic variance and at large due 
to pixel noise, with the fwhm chosen to be 10' {ig = 450) and a noise level 
of 6^K per 10' pixel (and with fsky = 1)- These error curves are not even 
visible in the £ range of the lower panel. These values are consistent with 
what might be expected from a very high precision satellite experiment like 
COBRAS/SAMBA |,|. 

What will limit this rosy picture is our ability to subtract foregrounds. Ul- 
timately, it will probably require a sophisticated combination of spectral and 
angular information, and cross-correlation with other datasets, such as X-ray 
and HI maps. With enough frequency bands covered, the prospects for sep- 
aration on the basis of spectrum alone is good. Figure || draws together the 
spectral signatures of the different sources of anisotropy that are likely to ap- 
pear and compares them with the frequencies that various experiments probe. 
Although the different signals are gratifyingly different, many parameters must 
be fit, either pixel- by-pixel, or using as well the different angular patterns that 
the signals will have. For example, extragalactic radio sources will have syn- 
chrotron spectra and a projected white-noise spectrum Cres.e like that 
shown in Fig. ^ for primeval galaxies; just as for the primeval galaxies, there 
could also be clustering contributions. The primeval galaxy frequency spec- 
trum would be similar to that of a cold dust component because of redshifting. 
The angular power spectra of Galactic bremsstrahlung and dust appear to 
obey Cres/ £~^, i-e-, fall to high i faster than scale invariance, apparently 
becoming small in the all important £ ^ 100 — 500 range, especially in the 
frequency range around 90 GHz [p^ [l^ . Complications will arise however, the 
most important being the non-Gaussian nature of the residuals and the multi- 
component nature of the dust, in particular the possible presence of cold dust 



2 CMB Distortions and Energetic Constraints 

In this section, I review the impact of spectrum observations on structure 
formation issues. We know from FIRAS that the CMB is well fit by a blackbody 
with T « 2.728 ± 0.004 K over the region from 5000Aim to 500/im [0, a 
number compatible with the 1990 COBRA rocket experiment of Gush et al. 
covering the same band, and also with ground based measurements at 
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Figure 2: The flat spectrum in thermodynamic temperature predicted for 
primary anisotropies is contrasted with the spectral signatures for other sources 
of anisotropy (normahzed at 4 mm): SZ anisotropies (long-dashed, with a 
sign change at 1300 /im); bremsstrahlung (short-dashed); synchrotron (dotted), 
with index vvarying from ps = 0.5 to 0.9 (intensity I^, oc i^~^^"); and dust (with 
index = 2 as indicated by a two-temperature fit to COBE), both the usual 
Galactic dust at 20 K (heavy solid) and dust at 6 K and 4 K (light solid 
lines, which could represent a cold Galactic component or, e.g., 30 K dust 
radiating at redshift ~ 5); a shallower, less physically-motivated, ad = 1.5 
dust opacity law for the 20K grains is also shown, appropriate for the single- 
temperature COBE fit. The frequency bands which various experiments probe 
are indicated, in particular the proposed bands for the NASA (MAP) and ESA 
(COBRAS/SAMBA) all-sky satellite experiments. There is a minimum of the 
Galactic foregrounds at about 90 GHz, the highest frequency COBE channel. 



centimetre wavelengths — although there is still room for significant spectral 
distortion longward of 1 cm. We must rely on indirect arguments based on 
primordial nucleosynthesis to constrain exactly when this photon entropy in 
our Hubble patch came into being, and whether this injection of energy would 
have a direct impact on short-distance structure formation. Energy injection 

/ s-0.39 

prior to zpi w 10^'^ I rish^/O.Ol ] is redistributed into a Planckian form: 

zpi defines the redshift of the cosmic photosphere. There have been heroic 
efforts to explain the CMB as starlight processed through exotic forms of dust 
(in particular long conducting needles) that would have happened much later. 
The observables from these models have never been fully worked out (see, e.g., 
but they are severely challenged by the absence of spectral distortions and 
the high degree of anisotropy required. 

Between zpi and zbe ~ 10'^-^(r2Bh^/0.0l) injected energy is redis- 
tributed into a Bose-Einstein shape characterized by a chemical potential, 
which FIRAS constrains to be [0 \n^\/T.y < 0.9 x IQ-^ (95% CL), trans- 
lating to a limit on energy of SEbe / Ecmb ^ 6.4 x 10~^ in this redshift range. 

Below Zy « 10^(r2Bh^/0.0l) the Compton ?/-distortion formula holds, giv- 
ing a unique signature to distortions, negative for frequencies below 218 GHz, 
positive above, and a stringent limit on the Compton-cooling energy loss from 
hot gas, SEcompton cool/ Ecmb = 4y < 6.0 X IQ-^ (95% CL). If there is no 
recombination, there is a constraint from the y-distortion on how early reheat- 
ing of the Universe can have occurred: Zmax,reh ~ 10'^ ''(r2Bh/0.02) '^^^Qni^ 
Ipl] , [2^ , but it is not very restrictive for the low Ob favoured by standard 
Big Bang nucleosynthesis and can be avoided if one can sustain a temperature 
of the cooling electrons to be nearly the CMB temperature. 

Compton cooling has been observed in more than two dozen massive clus- 
ters of galaxies above the 5-a level, including one at redshift 0.545, which tells 
us that the CMB existed by at least that redshift. With likely experimental 
sensitivity increases, the SZ effect may eventually offer a more powerful probe 
of the cluster distribution than X-ray observations do: instead of being a pro- 
jection of the square of the baryon density, it is a projection of the electron 
pressure and the decrease in signal with redshift is only a consequence of cluster 
evolution, not dimming by distance. Combining the SZ and X-ray observations 
is one of the main paths to Hq (and in principle qq), but is so far more confusing 
than enlightening, with values ranging from small {e.g., Hq = 38± 17 for Abell 
2218) to large (74 ± 29 for COMA); the hope is that a well-selected sample of 
clusters may help to reduce the biases. 

The integrated contribution of Compton-cooling from clusters and groups 
is not expected to be large for models of structure formation that reproduce 



the cluster X-ray temperature distribution function. For example |23 , for 



variants of adiabatic dark-matter dominated il = 1 models with ilp <C O, h 
0.5 and nearly scale-invariant initial conditions, the estimated value depends 
sensitively on erg, the linear amplitude of density fluctuations on the cluster- 



scale 8 h ^Mpc, and somewhat on the local curvature of the density fluctuation 
spectrum on cluster-scales: with as = 0.7, a hot/cold hybrid model with ili, = 
0.3 gives y ~ 0.3 X 10-^(17Be///0.05), a tilted CDM model with Hs = 0.8 
gives ~ 0.5 X 10~^(ilBe///0.05), with a similar value obtained for a = 1 
model. Lowering erg gives values well below 10^^, but also not enough high 
temperature clusters; raising it gives too many high Tx clusters (raising erg 
to unity still gives only y ~ 1.6 x 10^^(r2Be///0.05)). Here flseff takes into 
account the possible segregation of baryons from mass in clusters, modifying 
the value to be used over the primordial Qb- There is an even smaller effect 
associated with nonlinear Thompson scattering from the hot gas in the moving 
clusters. Nonetheless, the non-Gaussian pattern of Compton-cooling secondary 
anisotropics is accessible to experiment and will be a foreground to remove in 
future CMB anisotropy experiments 

The FIRAS limit on general secondary backgrounds (without a unique sig- 
nature hke BE or y distortions) is SE/E^rnbi^OO - SOOO^im) < 0.00025 {la 
CL). If pregalactic dust, or dust in primeval galaxies, exists, it will absorb 
higher frequency radiation (UV and optical) and down-shift it into the in- 
frared {e.g., ||2^ ); combined with the redshift, a sub-mm background is ex- 
pected but, with FIRAS, is now quite strongly constrained. The radiation 
could be largely shortward of 500/im; the peak in the vli, curve occurs at 
Apfc w 96(1 -I- z){30K/Td) for = 2 dust, which could be around 200/im if the 
dust is hot (which seems reasonable) or the redshift of bulge/elliptical forma- 
tion is low. The FIRAS constraint then applies only to the tail of emission. 
There is a tentative identification of a sub-mm background in the FIRAS data 
jlSf in the range 200 - 1000 ^m, with energy SE/Ecmb ^ 10"^ longward 
of 400 /im, which partly mimics the Galactic contribution (and could be 
partly due to cold high latitude Galactic dust ||l^). There are also residuals 
after source subtractions in the DIRBE data which could be interpreted as a 
cosmological infrared background at shorter (~ 1 — 200 fj,m) wavelengths at 
the 6E/Ecmb ^ 10^^ level These constraints on energy injection should 
be contrasted with plausible sources: For example, the nuclear energy output 
of stars with efficiency enuc radiating at redshift with an abundance fi* rel- 
ative to the CMB is E^/E^mb « 0.03(17*hV0.00l) [5/{l + z^)] [e„„c/0.004]. 
Massive stars in the 10-100 Mg range have e„„c ~ 0.004(M/100 M©)"-^, sat- 
urating at 0.004 above > 100 Mg (the Very Massive Object range). E^, has 
been used to constrain the role pregalactic black holes from VMO precursors 
could have played as dark matter. The massive star E^ is also tied to the 
heavy elements Z^jM they eject in supernova explosions. If the supernovae 
contribute a mean metal fraction Z to a gas of density flgas, Ep^eSN */ Ecmb ~ 
O.OOO8[Z/lO-3][^^^^^hVO.Ol][ZeJ/O.2](M/2OM0)°■^ Relaxation of the strin- 
gent energetic constraints is possible if either the energy was not reprocessed 
by dust (and so would reside in the near infrared where the DIRBE constraints 
are not nearly as strong j2^) or the dust was so hot that even with redshift 
effects it was shortward of 500/im. 

As Ikeuchi and Ostriker originally emphasized, a predominantly hydrody- 



namic explanation for cosmic structure development is a perfectly reasonable 
extrapolation of known behaviour in the interstellar medium to the pregalac- 
tic medium. However the Compton cooling limit constrains the combination 
fexp h^^Mpc) of filling factor f^xp and bubble formation scale Rexp to 

be < [(flBh/0-02)17y^^]-i {e.g., §). Further, if supernova explosions were re- 
sponsible for energy injection, one expects that the presupernova light radiated 
would be much in excess of the explosive energy (more than a hundred- fold), 
which would lead to much stronger restrictions on the model; and if the su- 
pernova debris is metal-enriched, the allowed amount of metals poses an even 
stronger constraint. (One may also argue that the specific tapestry that we see 
is too close to what straightforward gravitational instability predicts to war- 
rant consideration of a purely hydrodynamical model; i.e., that the explosive 
effects might be largely masked by subsequent gravitational instability. What 
does seem inevitable is that there will be a more limited local hydrodynamics 
role around collapsed objects.) 

3 Theoretical Issues and Sample Primary Power 
Spectra; COMBA 

The development of spectral distortions or angular anisotropics in the mi- 
crowave background is described by radiative transfer equations for the photon 
distribution function, which are coupled to Einstein's equations for the gravi- 
tational field and to the hydrodynamic and transport equations for the other 
types of matter present. The primary spectra are calculated by solving for each 
mode G {adiabatic scalar, isocurvature scalar, vector or tensor} the linearized 
Boltzmann transport equation for photons (including polarization) and rela- 
tivistic or light neutrinos, coupled to the equations of motion for baryons and 
cold dark matter, and the perturbed gravitational metric equations, possibly 
in the presence of vacuum energy or mean curvature. This is a well developed 
art the techniques for which have been described elsewhere {e.g., Q]) and will 
not be elaborated upon here. In homage to the high precision future that 
CMB experimentalists will provide for us, a large consortium of theorists who 
have developed computer codes to attack this problem fully or in various fast- 
computation approximations have combined under the acronym COMBA to 
deliver accurate validated calculations of Cis to the cosmological community 

I now sketch the Ci terrain in inflation-inspired models as the parame- 
ters defining the structure formation model are varied. Samples are shown 
in Figs. The "standard" scale invariant adiabatic CDM model (0 = 1, 
Us — 1, h — 0.5, fls = 0.05) with normal recombination shown in Fig. ^ and 
repeated in each of the panels of Fig. || illustrates the typical form: the Sachs- 
Wolfe effect arising from gravitational potential fluctuations dominating at low 
£, followed by rises and falls in the first and subsequent Doppler (or acoustic) 



peaks, arising from a combination of photon compression and rarefaction and 
electron flow at photon decoupling, with an overall decline due to destructive 
interference across the photon decoupling surface and damping by shear viscos- 
ity in the photon plus baryon fluid. A CDM model with very early reionization 
(at z > 200) shows no Doppler peaks, a result of destructive interference from 
forward and backward flows across the decoupling region, illustrating that the 
"short-wavelength" part of the density power spectrum can have a dramatic 
effect upon Cf, since it determines how copious UV production from early stars 
was. Lower redshifts of reionization still maintain a Doppler peak, but are 
suppressed relative to the standard CDM case (as illustrated by the Zr^h ~ 30 
model in Fig. |l] and in Fig. |^(e)). 

Figs. I^J^ include adiabatic scalar and tensor contributions. The relative 
magnitude of each is characterized by either the ratio of the quadrupole powers, 

Tts = d'p jcf^ , or the ratio of the dmr band-powers fts = {^f^)dm.rl (^i"^\mr- 
For the scale invariant cases, rts is taken to vanish. 

A simple variant of CDM-like models is to tilt the initial spectrum. The 
scalar tilt — Us — Xva defined in terms of the index , which is one for scale 
invariant adiabatic fluctuations. There is a corresponding tilt which charac- 
terizes the initial spectrum of gravitational waves which induce primary tensor 
anisotropics, vt = nt -f 3, where nt is —3 for a scale invariant spectrum. In- 
flation models give ft < and usually give < 0. For small tensor tilts, 
Tts ~ — 6.9i^t and vta ~ 1.3f(s are expected (with corrections given by eq. ^o|). 
For a reasonably large class of inflation models vt ~ Vs, but in some popular 
inflation models vt may be nearly zero even though Vs is not. Fig s. |l| and 1(a) 
show cf^ + derived for tilted cases when vt — i^s is assumed to hold. 

Fig. shows explicitly the contribution that C^^^ makes in one example. The 
tilt indices, especially I's, can also be complex functions of k in inflation models 
in which scale invariance is radically broken {e.g., |2^ ), in which case the Ce 
reflect the added complexity. 

Spectra for hot/cold hybrid models with a light massive neutrino look quite 
similar to those for CDM only, as Fig. ||(f) shows Q ||, |3|. This is true even 
for pure hot dark matter models [Q. 

The dotted Ci in Fig. ^ also has a flat initial spectrum, but has a large 
nonzero cosmological constant in order to have a high Hq, in better accord 
with most observational determinations. As one goes from £—2 to £ —3 and 
above there is first a drop in Ci |p8| , a consequence of the time dependence 
of the gravitational potential fluctuations (the integrated Sachs- Wolfe effect). 
Other nonzero A examples are given in Fig. ^(d). 

Open models like the Hq = 60 one shown in Fig. ^ have a nontrivial late- 
time integrated Sachs- Wolfe effect, like the f^A 7^ models do, but there is 

also a direct effect of the curvature on the mode function evolution, which 

1/2 

serves to focus the structure to a smaller angular scale (~ flnr ) than in the 
flcurv = 1 — r^tot = case. Of course whatever mechanism generated the 
ultra-large-scale mean curvature may well have had associated with it strong 



fluctuations on observable scales, so much so that this is an argument against 
large mean curvature because of the absence of such effects in the CMB. Even 
if the background curvature is determined by an entirely different mechanism, 
it should influence the fluctuation generation mechanism. An open issue in 
open models has always been what is a natural shape for the spectrum for k 
near d~^^^ = Hq^IL^v Power laws in kdcurv, \J {kdcurv)^ — 1 etc. have often 
been adopted but if the fluctuation generation mechanism is quantum noise 
in inflation, there is a natural adiabatic spectrum expected which is a simple 
generalization of the nearly scale invariant spectra of ftcurv — inflation models 
p9| , |30| , Fortunately for £ ^ 2fl^^rtlLrv this issue is not a factor, so that 
intermediate and small angle predictions are relatively unambiguous. 

Inflation-based models with isocurvature rather than adiabatic initial condi- 
tions are strongly ruled out by the CMB data if they are nearly scale invariant 
p^ , |3^ , but could contribute at a subdominant level to the adiabatic fluc- 
tuations. Even allowing for arbitrarily broken scale invariance in the initial 
fluctuation spectra, the allowed region for pure isocurvature baryon or CDM 
models has been shrinking fast as the data has improved. 

Defect models have (knot-like or string-like) localized topological field con- 
figurations acting as isocurvature seed perturbations to drive the growth of 
fluctuations in the total mass density. On large angular scales the defect models 
lead to a similar nearly scale-invariant spectrum for Ce as for inflation-inspired 
adiabatic perturbation models |34, 36, |8[ On smaller scales, the 



spectra are sufficiently different from adiabatic inflation-inspired spectra to 
sharply test these competing pictures of cosmic structure formation in the next 
generation of CMB experiments. The non-Gaussian nature of defect-induced 
anisotropics also adds another point of differentiation among the models. 

The angular power spectra generally differ substantially as a function of 
multipole £ so that even with the experimental results we expect in the very 
near future major swaths of cosmological parameter space can be ruled out. 
However, different combinations of the parameters n^, r^s, ^^A etc. can lead 
to nearly identical spectra [D, as Fig. ||(f) illustrates. Superposed upon the 
spectra in Fig. ^ are theoretical band-powers derived for a variety of anisotropy 
experiments. Fig. || also shows 10% 1-sigma error bars: COBE achieved 14% 
errors with 4 years of data; to achieve this with smaller angle experiments one 
needs to have about the same number of pixels as COBE, but scaled to the 
beam size hence covering a smaller region of the Universe. So far none of the 
smaller angle data sets have the 650 or so fwhm-sized pixels COBE does, but 
this is the stage we are now entering [|9j. 

As noted before, even if there were idealized perfect all-sky coverage with 
noise- free versions of the experiments of Fig. there would still be cosmic vari- 
ance errors on the band-powers, but these are ^ {£) ^ pO| , much smaller than 
the size of the points. Thus it appears that by using (perfect) CMB experiments 
which are sensitive to a wide range of angular scales, we can distinguish even 
among the nearly degenerate theoretical models shown, and be able to measure 



the parameters that define the variations in these models. A closeup view of 
examples of how we can measure very fine differences in models is shown in 
Fig. ^, using the detector sensitivities and long observing times that satellite 
experiments now currently feasible can achieve M, @|. 



4 Secondary Anisotropy Sources 

Reionization and Primary CMB Anisotropics: An important issue as- 
sociated with the early energy injection described in § ^ is what impact it will 
have on the primary anisotropics of the CMB. The physical processes impor- 
tant in the recombination of the primeval plasma have been well understood 
since shortly after the discovery of the CMB. The comoving width of the re- 
gion over which decoupling takes place if there is normal recombination is only 
(5 — 10)r2„r h^^Mpc, where Qnr is the density in non-relativistic particles 
(CDM, baryons); the viscous damping scale of the photon-baryon fiuid prior to 
recombination is slightly less. The associated angular scale can be character- 
ized by a multipole number, £ ~ 1000, above which anisotropics are strongly 
damped (Fig. [l]). The associated natural 'coherence' angle, ~ 10', defines which 
experiments are most useful to do if we wish to probe the moment when the 
photons were first released to freely propagate from their point of origin to us, 
without much further modification, apart from some gravitational redshifts, 
some lensing, and possibly some scattering from hot gas. 

The main effect that reionization of the Universe has on anisotropics is 
to lower their amplitude by a factor AT/T oc e~^^, where is the optical 
depth to Thompson scattering. If Zreh is the reionization redshift and Zi^^^^i ~ 

102 i(l7sh/0.02)"^/^ nl/r^ is the redshift one would need to reionize by to get 
a Thompson depth of unity, then Cc = [(1 + 2)/(l + With standard 

Big Bang nucleosynthesis values for fi^, getting Zfc=i much below 100 seems 
unlikely. Zreh is presumably the redshift by which the first nonlinear objects 
form in sufficient abundance to allow enough massive star formation to occur 
to cause pregalactic HII regions to overlap, a quantity largely determined by 
the short-distance density fluctuation power in the structure formation theory 
in question, but subject to many uncertainties: the entities which form may 
well be rather fragile with a small binding energy, easily disrupted by the 
massive stars they generate; on the other hand, the amount of nonlinear gas 
could be amplified by the explosion of such stars sweeping up shells of gas 
far from the parent object. Thus z^e/i depends upon how rare the ionization- 
generators can be. For infiation-based CDM models and variants with light 
massive neutrinos or nonzero A, z^e/t values ranging from 5 to 60 seem plausible, 
hence Cc ranges from negligible to substantial, < 1/3, but not so large as to 
fully erase anisotropics, which would occur on scales below Tdec/T^, where Tdec 
is the horizon scale at photon decoupling, corresponding to ^ '--^ 100 for the 
"standard" CDM model. See Figs. |^ and ^(e) for Zreh 30, 50 examples. 
In isocurvature baryon models with (nearly) white noise initial conditions 
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Figure 3: (previous page) Spectra for a variety of inflation-inspired models, 
normalized to the COBE band-power. Theoretical band-powers for various 
experimental configurations are placed at horizontal error bars extend 

to the e~^^^Wmax points. Unless otherwise indicated, riBh^=0.0125, h=0.5, 
Us—l; when the gravity wave contribution is nonzero, vt—i^s is assumed (and 
Tts = /C2 is ^ —7vt). The untiltcd 7is=l, fts=^ model is repeated in each 
panel (solid line), (a) CDM models with variable tilt Ug. (b) Us — l models with 
r^sh^ changed, h fixed, (c) Us—l models, with fish^ changed, fis fixed, (d) 
ns—\ models with fixed age, 13 Gyr, but variable Hq and f2A = 1 — ^cdm — ^b 
(.92, .79, .43,0 for 100,80,60,50). (e) CDM models with very early reionization 
at Zreh ?J 150 (equivalent to 'no recombination'), and later reionization at 
■2re?i=30,50 are contrasted with standard recombination (SR). The Zre/i=50 
spectrum is close to the rts=0.95 spectrum with SR (thin, dot-dashed): the 
moderate suppression if 20 < Zreh i$ 150 can be partially mimicked by decreas- 
ing ris or increasing h. (f) Sample cosmologies with nearly degenerate spectra 
and band-powers. Dashed curve: increasing SIa is compensated by increasing 
h. Dot-dashed curve: tilting to 71^=0.94 (fts=0.42) is compensated by increas- 
ing ilA to 0.6. The dotted hot/cold model curves (with Vl^ indicated) are 
nearly identical to the standard CDM one, but even these few percent differ- 
ences can be distinguished in principle by satellite all-sky experiments with 
currently available detector technology. 



1) beam = 20', 20' pix; noise per pix = 

2) beam=10', 10' pix; noise per pix = 6/i,K 
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Figure 4: This shows the ability of satellites to measure cosmic parameters to 
high accuracy. The relative difference of the power spectrum in question from 
a comparison spectrum (both normalized to the 4-year dmr (53-|-90-|-31)A-|-B 
COBE maps) are shown so that the few percent deviations can be clearly seen 
over the entire i range. The lighter lines are 1 — sigma error bars for all- 
sky coverage (averaged over the smoothing width shown, with ^Aln^ = 0.05) 
and include cosmic variance (dominant at low £) and pixel noise at 20/iK or 
6/xK (dominant at high £), with the very rapid growth relative to the theory 
curve at high £ coming from the finite beam-size (with the fwhm indicated, 
corresponding to a Gaussian filter in multipole space of is = 404 and £s = 809 
respectively). The first choice corresponds to the NASA satellite experiment 
MAP, the second choice to the ESA mission COBRAS/SAMBA, assuming the 

— 1/2 

entire sky is usable (errors scaling oc 

fsky )■ The ultimate accuracy achievable 
will depend upon the decontamination of the primary signal of non-Gaussian 
Galactic synchrotron, bremsstrahlung and dust signals. The models shown all 
have a uniform age of 13 Gyr, Qcdm + ^mi^ + Q\ + = 1, f^sh^ = 0.0125, 
= 1 and no gravity wave contribution. Notice the scale change for the 
hot /cold model panels. (One species of massive neutrino was adopted for these 
two cases.) 



popular in the late seventies, the first objects collapse at z ^ 300, making reion- 
ization easy, and, indeed, expected. Thus, large viscous damping is expected. 
However, one can still have a peak in at ^ ~ 200 for open models, essentially 
because the modified angle- distance relation in curved universes shifts power to 
higher £. Early ionization seems plausible, but by no means certain, in models 
in which there are isocurvature seeds, such as in texture models [ p6[ . 
Reionization and Quadratic Nonlinearities in Thomson Scattering: 
These can sometimes dominate over the first-order anisotropics if the latter are 
strongly damped and there is early ionization. Even if there is early reionization 
in nearly scale invariant models, there is generally not sufficient power on small 
length scales for this Vishniac effect |Q to be important. Thus it can usually be 
ignored in inflation-based models. This is not so for isocurvature baryon models 
|]33| in which the initial spectral index n^s, considered a free parameter, is 
between -1 and on phenomenological grounds. Such a steeply rising spectrum 
implies short-distance effects are very important, and give predicted sizable 
signals in e.g., the ovro and VLA window of Fig. [l[ 

The Rees-Sciama Effect: In flat flnr — 1 models, the gravitational potential 
fluctuations have constant amplitude in the linear phase of evolution. Weak 
(or strong) nonlinear evolution induces time dependence in $Ar which induces 
anisotropy by an integrated Sachs- Wolfe effect. For inflation- inspired models 
this turns out to be a very small correction and can be largely ignored. 
The Influence of Weak Gravitational Lensing on the CMB: Another 
nonlinear effect (on the distribution function) is gravitational lensing which 
bends, focusses and defocusses the CMB photons as they propagate from de- 
coupling through the clumpy medium to us. Given the difficulties that as- 
tronomers have had detecting lensing, with the best observations coming from 
clusters of galaxies, it may seem obvious that the effect on the ~ 10' coherence 
scale typical for primary CMB anisotropies is likely be quite small; and this is 
what (most of) the people who have investigated the effect have found. Lens- 
ing conserves the total angular power, it just rearranges it, by smoothing the 
Doppler peaks. The typical range in £ over which the power is spread in Ai/£ is 
basically the weak-lensing shear, about 10% to 20% or so at a few arcminutes, 
depending upon the model |^^ ; this is in agreement with the levels estimated 
by people advocating using the influence of weak-lensing on the ellipticities of 
faint galaxy images to determine the mass density power spectrum e.g., \ ^7\ . 
Sunyaev-Zeldovich Fluctuations and the Moving Gluster Effect: The 
SZ effect has been discussed in § ||. The moving cluster effect is the nonlinear 
Thomson scattering of the CMB photons from plasma conflned to clusters, 
moving with it. It is predicted to be quite a bit smaller than the SZ anisotropy 
level, and as can be seen from Fig. |^ the contribution of SZ from clusters is 
small relative to the primary anisotropy power ||2^ . However, the distribution 
is non-Gaussian, concentrated in the pressure-peaks in the medium, especially 
the clusters. This means a search for the ambient SZ-effect (where the SZ 
sources are not known beforehand) could be very promising . 
Starbursting Primeval Galaxies: As discussed in § || on CMB distortions 



and energetic constraints, the prospects for seeing these with sub- mm telescopes 
probing tens of arcseconds (e.g. SCUBA filter in Fig. |l|) seem quite good 



5 Anisotropy Experiments and their Band-Powers 

The importance of the large-angle COBE dmr ^ detection for testing 
theories of cosmic structure formation can hardly be overstated. Smaller angle 
experiments are now also achieving results that can be combined with COBE 
to constraint model parameters. In this section, I review the experiments. 

The COBE data stream gives temperature differences in pixel pairs sep- 
arated by 60°, but an inversion yields 6 maps, for the A and B channels of 
the 3 frequencies shown in Fig. ||, 31, 53 and 90 GHz. Figure || compares an 
optimally filtered map of COBE's data with what a scale invariant = 1 dark 
matter dominated model would look like under the same filtering. 

The i = 2 power uses the 4-year quadrupole value [^9[ |5^ , determined from 
high Galactic latitude data. It is the multipole most likely to have a residual 
Galactic signal contaminating it, possibly destructively, and the "systematic" 
error, the dashed addition to the statistical error bar (solid), reflects this. The 
two heavy points at ^ ^ 7 are band-powers derived for the 4-year dmr 53-f 90-1-31 
GHz maps |l|, |l^ , the solid point assuming a z^at = spectrum, 

the open marginalizing over all possible vat, where v^t is defined by eq.(|^). 

The two points at ^ ~ 10 arc for the firs map ||5^, solid with the restriction 
z/AT — 0, open with i/^T allowed to float. The coverage of the firs experiment is 
much less extensive and more inhomogeneous over the observed patch than for 
dmr. it was a balloon experiment (taking useful data for only about 5 hours) 
with bolometer detectors probing the 4 frequency channels shown in Fig. ^ 
Only the 170 GHz channel has been fully analyzed. The pixel size, 1.3°, and 
the beam-size, is ~ 34 (a 3.9° fwhm beam) are half the COBE values. The firs 
map has a positive cross-correlation with dmr . The band-powers of dmr and 
firs are quite comparable, almost independent of degree of the signal-to-noise 
filtering and which frequencies are probed. And it is only weakly dependent 
upon the slope vat (as the formula for (C^}^;,^^ as a function of vat given by 
eq.(|l2|) shows). 

The effective slope of the standard CDM model is vat ~ 0.15 over the 
dmr band; variation in fi^ and Ho does not change this very much, dmr and 
firs have enough coverage in £-space that one can estimate the spectral index 
from the data as well by Bayesian means. If there is no filtering the 53+90+31 
GHz 'A+B' index is vat is 0.07 ± 0.28 for dmr, firs gives a steeper index, 
O.B^Qg, but the one sigma error bar encompasses and there is clearly a small 
angle residual 'noise' driving the higher values |l|. When CDM models are 
marginalized over erg the preferred dmr index is rig = 1.02i:g for the case with 
no gravity waves {vt — 0), and Us = 1.02^'^g with gravity waves (vt = Vg) 
|ll3|. Band-powers for specific £ ranges also show the nearly flat character for 



Ce (light open points at £ ~ 4, 8, 16 from |51 ). 



4-yr COBE North Galactic Pole Map 

UNFILTERED OPTIMAL FILTER 




140 degrees ^ ^ 140 degrees 

(b,c,d) contours: AT/T = ±0.55x10-^ n, n^l 
(a) contours: AT/T = ±0.55x10-5 2"-i, n^l 

Figure 5: 140° diameter maps centered on the North Galactic Pole are shown 
for a realization of a CDM C^-spectrum convolved with the dmr beam in (c). 
No noise has been added. This is how the primary sky would appear in a 
Us = 1 CDM Universe with cg = 1.2 (or in a flmi^ = 0.2 hot/cold universe with 
(Ts = 0.8), the most likely values for the dmr data. This is contrasted with the 
4-year dmr (53+90+31)a+b data map shown in (a) and the map after the data 
has undergone optimal signal-to-noise filtering in (b) (using the same C^-shape 
and amplitude for the filter) . The statistically significant features are also seen 
in each of the dmr channel maps after optimal filtering (which preferentially 
removes high angular frequencies, more so for noisier maps). Thus, to compare, 
(d) shows the theoretical realization after passing (c) through the same optimal 
filter used for (b); the average, dipole and quadrupole of the full \b\ > 20° sky 
were also removed, an effective low ^-filter - if they stay in, the maps look 
similar to the unfiltered theory maps except small scale smoothing leads to 
loss of the higher contour levels. Note that the contours are linearly spaced at 
±15n fiK for all but (a), for which the spacings are ±15,30,60, 120 /iK. The 
maps have been smoothed by an additional 1.66° Gaussian filter. 



The Tenerife point at ^ ~ 20 uses combined 15 and 33 GHz data. The 
amphtude is compatible with dmr and firs, with CDM-hke models, and also 
have features in common with dmr. 

We now come to the crowded region from two degrees to half a degree. The 
next two experiments are the South Pole HEMT experiments of the UCSB 
group and the Saskatchewan (Big Plate) HEMT experiment of the Princeton 
group. The lower open circle is from a joint 4-channel analysis of the 9 and 
13 point sp91 scans |5^, ^ |l3j (with the individual 9 point and 13 point 



OCT solid point is for a simultaneous 
5q, Esf, with separate values for the 



values given in the lower panel). The upi 
analysis of all channels of the sp94 data 
Ka (~ 30 GHz) and Q (- 40 GHz) HEMT bands (Fig. g) in the lower panel. 
The solid triangle in the upper panel is the sk93 result the big solid 

circle at ~ 80 in the lower panel is the sk93+94 result (with calibration 
uncertainties adding another 14% error to the statistical error shown). The 
nearness of the sp94, sk93 and sk93+94 band-powers, and the demonstration 
for both experiments that the preferred frequency dependence is nearly flat in 
AT/T and many sigma away from bremsstrahlung or synchrotron, the expected 
contaminants in this 30-40 GHz range, lend confidence that the spectrum in 
the £ ^ 60-80 region has really been determined; and it looks quite compatible 
with the COBE-normalized CDM spectrum: sp94 gives as = 1.261;^^, and 
sk93+94 gives as = 1-211^11 lH' very close to the dmr value as = 1.20 ± 0.08 
and the firs value erg = 1.27 ± 0.30. The 5 heavy open circle points probing ^'s 
ranging from 60 to 400 repeated in the upper and lower panels labelled sk95 
are combined sk94+95 results |5|]. The estimated 14% error in the overall 
amplitude because of calibration uncertainties are included. The large ^-space 
coverage from this one intermediate angle experiment gives a first glimpse of 
the £-space coverage that will become standard in the next round of anisotropy 
experiments. 

Python |6^, py, the heavy solid curve at ^ 90, is sensitive to a wide 
coverage in ^-space as the horizontal error bars in the top panel indicate. Argo 
jsH , ar, a balloon-borne experiment, is next. The next five points in the 
lower panel are from the fourth and fifth flights, M4,M5, of the MAX balloon 



experiment 1 63 , 6^ . Because the filters changed with frequency, the points are 
placed at the average over all max filters. In the upper panel three max4 scans 
are combined into one data point as are two max5 scans. The lines ending in 
triangles at ^ « 145 and 240 denote the 90% limits for the MSAM [|4| single 
{msam2) and double (msamS) difference configurations. A limitation on these 
balloon experiments is the ^ 5 hours over which data can be effectively taken. 
Planned long duration balloon flights that would circle Antarctica for about a 
week would allow extensive mapping at high precision to be done, and a number 
of groups have been proposing designs {e.g., ACE, Boomerang and Top Hat). 

The CAT points at £ ^ 400 and 600 represent a very different ex- 
perimental technique, interferometry. CAT is a 3-element synthesis telescope, 
probing ~ 15 GHz frequencies with a 27' synthesized beam and a 2° field- 
of-view (the fwhm of the individual telescopes). It is a precursor to the larger 



VS A (Very Small Array) , covering a wider frequency range with more telescopes 
and a larger (4°) FOV. Two other CMB interferometers are also planned: CBI 
and VGA. The ovro experiments also probe radio frequencies, but using single 
dishes. The historically important 1987 ovro 7 point upper limit shown 
used a 40 meter dish. Detections using as well a 5 meter dish have now been 
found with ovro and give a value in between the 2 CAT points with about the 
same amplitude. The open triangle at ^ « 160 denoting the 95% credible limit 
for the sp89 9 point scan was also historically important. WhiteDish 

pSf had a small amplitude filter, a hint of a detection in the m = 1 mode, and 
a 95% limit in m = 2 mode at ^ w 520, wd2. 



6 COBE-normalization of Post-Inflation Fluc- 
tuations 

For early universe calculations and also to characterize the initial conditions 
for the photon transport through decoupling, the power in adiabatic scalar 
fluctuations on scales beyond the Hubble radius is best characterized in terms 
of quantities which become time-independent. Some examples are the spa- 
tial curvature of time surfaces on which there is no net flow of momentum 
(oc (ficom), as Chibisov and Mukhanov emphasized long ago |^9[ ^ |7l), the 
expansion factor fluctuation, Slnaln*, on time surfaces with uniform space 
creation rate H^, ^ and Bardeen, Steinhardt and Turner's C ifl, 13- 
An initially scale invariant adiabatic spectrum has fc-independent power per 
d\nk in these variables (for k/{Ha) ^ 1), while for models with spectral tilt 
Vg, we have Vip^^^^k) — 'Pip^^^{Tf^^){kToY% where we use the instantaneous 
comoving horizon size at the current epoch, tq, as the normalization point. For 
CDM-like models (those with fl — Qnr — 1 and tq = 2Hq^), these are related 
to the portion of the dmr band power (Ci)^^^ in the scalar adiabatic mode, 

fdmr ^ (^e)dmr/ (l+'^ts); and to the quadrupole power, 

by 

'P^..,Aro^) ~ 23.5(cf g-i.99..(i+o.i..) ^ 23.6Cf' e-'-'^' , (8) 

i.e., about 3 x 10^^. This relation is very insensitive to variations in h and 
ri_B. For scales of order our present Hubble size, we also have Vc; ~ ir^*jv ~ 
^'P(Sp)hoT^ where $Ar is the perturbed Newtonian gravitational potential and 
{5p)hor is the density fluctuation at 'horizon crossing', defined by fcr = 1. 

Quantum noise in the transverse traceless modes of the perturbed metric 
tensor would also have arisen in the inflation epoch and for many models may 



have been quite significant |73, [7q, 77 . The gravitational radiation power 



spectrum Vqw = T^h+ + 'Phx is the sum of the two independent gravitational 
wave polarizations. It is related to the amplitude of the dmr band power 



('^e^^)drnr = C^^) dmr'^ts / + fts) and to the quadrupole C^^^ by 

VcwiTo') « 17.6(cf \_ e~'-92Mi+o.i..) « i3.4cf ) e-i-^^-^* , 
T'GwCfc) = ^Gw(ro-')(A:To)''* . (9) 

The inflation model determines the ratio of Pgw{tq^) to 'Pv'com ("^(T^)- I* is 
generally related to the tilt of the gravity wave spectrum, and this can in turn 
by used to relate the ratio of dmr band-powers (and quadrupoles) to the tilts 

(for r^^fjc — ^curv — 0)- 

Vgw{to^)/V^^^^{t^^) = -^vtl{\ - vtl2) ; hence 

\W Idmr ^ ' ' 

*^-cf) (l-^*/2) ■ ^ ^ 

The tensor tilt is simply related to the deceleration parameter q = —ad /a?' of 
the Universe in the inflationary epoch, vt/2 Ki 1 + q~^-^ although 1 + q~'^ is the 
leading term for the scalar tilt, other terms can dominate when the deceleration 
is near the critical deSitter-space value of —1 {e.g., |7^). Thus although vt is 
negative, Vg may not be. 

When assessing the effect of gravity waves on the normalization of the 
spectrum, it is useful to consider two limiting cases: Vg ~ t^t, which holds 
for the widest class of models, including power law and chaotic inflation, and 
Vt ~ 0, with Vg arbitrary, which holds for some models such as "natural" 
inflation. 

There are also corrections as one goes away from the D,nr — ^tot — 1 
models. For example, models with nonzero cosmological constant Qvac, but 
rinr + ^vac = 1, have Pqw / {Cg /^^r being only weakly dependent upon fl^ac 
whereas T^c/ i'^i^^) dmr °^ (1 ~ ^-^^lac) strongly dependent upon it. 



7 Using COBE, FIRS, SP94, SK94, ... to fix as 

Before the COBE detection, normalization of the density spectrum was done 
using (78, the rms (linear) mass density fluctuations on the scale of 8 h~^Mpc, or 
to a biasing factor bg for galaxies, which was usually asstimed to obey bgas « 1 
e.g., j?!]. The relation of as to the initial power spectrum amplitude 
'Pip^^^{tq ) is more sensitive to the specifics of the model, such as type of 
dark matter, spectral slope, h, than is the dmr band-power relation given in 
§ 1^. Comparing erg estimates from large scale structure observations with the 
COBE-normalizcd value is thus extremely important for constraining cosmo- 
logical parameter space. In this section I present a useful functional form. 



ag,{T,iys,vt,rtcdm,^m}^,^vac), where F parameterizes the density power spec- 
trum, then use constraints on erg, F and the tilts to sketch which models can 
already be ruled out. 

A byproduct of the linear perturbation calculations used to compute AT/T 
is the transfer function, which maps the initial density fluctuation spectrum 
in the very early universe into the final post-recombination one. Many fits to 
transfer functions have been given in the literature. One of the most useful 
exploits the approximate scaling with the "horizon" at redshift flnr/^er when 
the density in nonrelativistic matter, flnrO,^^ , equals that in relativistic matter, 
^7era-^ fc^eq = 5 F"^! h-^Mpc. The factor Fe, = r2„r h [17^ J(1.6817^)]~i/2 
provides the main shape dependence, but a further correction factor can ap- 
proximately incorporate the effect of baryons for Qb ^nr over the large scale 
structure region in fc-space ^ . One functional form for this is 

Vp{k) oc fc4+.= W |i ^ [ak+(hkf'^ + (ckfYy^'" Dl , (11) 



(a, 6, c) = (6.4, 3.0, 1.7)T-^ h^^Mpc ,p ^ 1.13 
ft 



(1.68f7.^) 



g-(OB(l+J2-^i(2h)i/2))-0.06) 



(Another functional form uses the well known fls BBKS version of the 
CDM transfer function as a base upon which the F variations are imposed.) 
Generally, more scales are needed to characterize the spectrum than just kfjeq] 
e.g., the collisionless damping scale for hot dark matter (massive neutrinos) 
KLmp ~ 6 (f7„^f]42h)2)-i/2(^^^/2)i/2 h-iMpc, with g™^ the number of 
massive species (counting particle and antiparticle) . The correction factor for 
massive neutrinos, -D^, is fit by js^, and is quite accurate even for finite Q,b 

il- 

Estimations of ag from the dmr data for selected models can calibrate a 

1 /2 

scaling relation between erg and {Cf)J^^ found using the "naive" Sachs Wolfe 
formula AT/T ^ <l'jv/3 relating temperature fluctuations to gravitational po- 



tential fluctuations |82H5|, [78[ |2| 



1.25 10^(C,)1^, »-r"(2(r - 0.03)) ,.63.. 
r-law. - 7^ (1 + (1 + o.55(f7_/0.3) V2) ^ ' ^^2) 



4yr(53&90&31)a&b: lO^{Ce)li^.^ ~ [0-82 + 0.26(1 - '^'f'^] x 11;°^ 

fsw « (1 + 0.12r!„,e)(l + ^l°ac) , ''AT ~ 0.15(1 - n,ac) + "s , 
(L ~ Vtl I) 

In Figure ^, the top left panel shows the average and ±1ct variation of as 
against tilt for a pure CDM model when no gravity wave induced anisotropics 
are included (upper hatched region) and when they are, for equal tensor and 
scalar tilts (lower hatched region). The formula labelled H + C shows how 



-.06 ' 



much (Tg is reduced for hot/cold hybrid models (for one species of massive 
neutrino, the other two assumed to be effectively massless). The heavy closed 
circles with the error bars denote values obtained in a Bayesian analysis of the 
4-year dmr (53+90+31)A+B COBE maps using the exact Ci. The upper right 
panel shows the variation for hot/cold hybrid models with variable f2„ii/, for 
the untilted (upper) case and the ?is=0.85 tilted (lower) case (with no gravity 
wave induced anisotropies). The open circles with error bars to the left of the 
dmr closed circles denote erg's derived from the sp94 data, while the squares 
denote erg's from the sk93+94 data. 

The next four panels show erg against the Hubble parameter h for fltot = 
1 models defining a fixed timeline in parameter space; i.e., enough vacuum 
energy, ^l^ac = A/(87rGAr) = 1 — has been added to keep the total unity 
and the age constant. Using the iluac dependences of fswi ^at and fts allow 
good (Tg fits.0 All models have no mean curvature, 17Bh^=0.0125, with the rest 
of the nr-matter in cold dark matter {^nr = ^cdm + ^b)- For these sequences 
of models with a uniform age tg, the variation of fl^ac with Hubble parameter 
(the rising curve in Fig. ^ is 



h = hi^,J 1 , /ii = 0.5(13Gyr/io) 



n^acih) ~ 0.9(0.3(h/hi - l)"-3 + 0.7(h/hi - 1)0-4) (-13) 

The latter is a rough inversion. The ages shown in Fig. ^ bracket a recent 
Hipparcos-modified estimate for globular cluster ages, 11.7 ± 1.4 Gyr |8^, with 
perhaps another Gyr to be added associated with the delay in globular cluster 
formation. The fl^ac = model with 13 Gyr age is therefore the Hq = 50 
standard CDM model; for the 15 Gyr age, Hq = 43 and for the 11 Gyr age, 
Ho = 59. 

The points with error bars again denote the dmr, sp94 and sk93+94 values. 
The 4-year dmr data alone places no useful constraint on Clyac or Hq: for a 
scale invariant spectrum with 13 Gyr age, ^vac < 0.85 and Hq < 88 at the Icr 
level. The constraint is slightly tighter but still weak when sp94 and sk93+94 
are added: flvac < 0.56 and Hq < 65 at Icr, < 0.83 and < 85 at 2cr. For a 15 
Gyr age, fi^ac < 0.35 and Hq < 51 at la, < 0.83 and < 73 at 2a. When the 
slope is marginalized over, the constraints relax. In the figure, the untilted 13 
Gyr case is augmented by a tilted case, with ris = 0.9 and a matching tensor 
tilt, = = —0.1. It is not quite flat enough for there to be problems 
matching both the sk and dmr data within the error bars {e.g., for Hq = 50 
and 13 Gyr, we get ||l3| = 0.95lo o8 with gravity waves, 0.93to i2 without; 
for Hq = 70, Us = 0.92l:°;o5 with, 0.88t°;i2 without; when Hq is marginalized. 

Although fsw = 1 takes into account some of the enhancements over the naive Sachs— 

1/2 

dmr 
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Wolfe formula by normalizing to the calculated ag-(Ci}, relation for standard CDM, it 



1/2 

does not take into account the enhancement of ) associated with the time dependence 
of the gravitational potential when A dominates: to take this into account, we need fsw to 
exceed unity, as shown. 



these numbers are O.OS^q'q^ and 0.92t°;i^ .) As is visually evident from Fig. hi, 
using the sk95 rather than the sk94 data with dmr^ increases the preferred 
Us- However, when all smaller angle data is used, and Hq is marginalized. 
Us = 1.0 ± 0.04 is obtained, a rather encouraging result for inflation-based 
models. 

The mass enclosed within 8h~^Mpc is that of a typical rich cluster, 1.2 x 
lO^^ilnr (2h)^^ Mq. Because rich clusters are rare events in the medium, their 
number density is extremely sensitive to the value of erg; the relative numbers 
of rich and poor clusters also depends upon the shape of Vp{k) in the cluster 
band, ~ (0.15 — 0.25) hMpc^^, i.e., on iis and T. Cluster X-ray data implies 
0.6 ^ (Ts < 0.8 for CDM-like Vlnr = 1 theories, with the best value depending 
upon r. Us, some issues of theoretical calibration of models, and especially 
which region of the drid/dTx data one wishes to fit, since the data prefer a 
local spectral index d\nVp/d\nk substantially flatter over the cluster region 
than the standard CDM model gives I believe a good target number is 

0.7 and values below about 0.5 are unacceptable, but because CDM spectra 
do not fit the data well, this normalization depends upon whether one focusses 
on the high or low temperature end. Other authors who concentrated on the 
low to median region found lower values for tIs = 1 models, 0.57 ± 0.05 js^ 
and 0.50 ± 0.04 but do not fit the high Tx end weU.0 In the top left and 



right panels, the two vertical lines denote two estimates of trg from clusters for 
^vac=0 (which depends somewhat upon tilt). For f^A > models, higher values 



are needed to fit the cluster data |82[ adopt ^nr '^^ the correction, 

psf give a more moderate dependence, ri~°-^'^+°-^'^^"'\ The rising curves with 
error bars in Fig. ^ using the ^nr '^^ scaling, show the higher and lower ag 
estimates. Allowed models would have to lie in the overlap region between the 
cluster (Tg and the dmr as . 

The upper rising regions also roughly denote the dg behaviour as derived 
from optical galaxy samples, in units of bgcrg where bg is the biasing factor for 
galaxies.^ How seriously we take the constraints derived using flows and redshift 
surveys depends upon how reliable we think the indicators are - a subject of 
much debate. Much work has also gone into relating bulk flows of galaxies to 
CMB anisotropy observations on intermediate angular scales, since they probe 

small upward correction should be applied to these low estimates to account for the 
nonzero redshift of the calibrating samples. 

*There are many estimates of the combination crgfij^,?'^ that are obtained by relating the 
galaxy flow fleld to the galaxy density field inferred from redshift surveys, which all take 
the form agfl^^^ = [bg(7g]/3g, where /3g is a numerical_factor whose value depends upon 
data set and analysis procedure. Strauss and Willick |p9|] have reviewed the rather varied 
estimations and give raw averages: 0.78 it 0.33 for IRAS-selected galaxy surveys, 0.71 it 0.25 
for optically-selected galaxy surveys. Even more recent estimates give /3g ~ 0.4-0.6. These 
as estimates rely on the simplifying assumption of a linear amplification bias bg for galaxies. 
The traditional estimate is bg = ^ , but bg can depend upon the galaxy types being probed, 
upon scale, and could be bigger or smaller than cr~^, and certainly cannot be determined by 
theory alone. A recent estimate using the Mark III velnrit y d ata is agD,^^^ ^ 0.85 it 0.1, 
with sampling errors adding another ~ 0.1 uncertainty M, blj. (The f^J^V?® is the factor by 
which the linear growth rate D/D differs from the Hubble expansion rate a/a.) 



the same bands in fc-space, and this has led to significant constraints, but is 
also subject to much debate. 

The shape of the density power spectrum is almost as powerful a constraint 
as as is. To fit the large scale galaxy clustering data, in particular the APM 
angular correlation function, requires 0.15 < r+z/s/2 < 0.3 |8^, assuming 
the galaxy power spectrum on large scales is a linear amplification of the (linear) 
density power spectrum. In the lower panels of Fig. ^, the dropping curve 
shows ri„rh and the (almost indistinguishable) dashed one shows F; the target 
range is shown by the line at 0.2 with a few error bars. To lower F + v^j^ 
into the 0.15 to 0.3 range one can ^ tilt the spectrum (expected 
at some level in inflation models); lower h; lower r2„r; or raise f2er (= LBSfi-^ 
with the canonical three massless neutrino species present). Raising also 
helps. Low density CDM models in a spatially flat universe (i.e. with A > 0) 
lower r2„r to 1 — JIa- CDM models with decaying neutrinos raise fier p3| : 
F w 1.0817„rh(l + 0.96(TOi/rd/kevyr)^/'^)~^/^, where is the neutrino mass 
and is its lifetime. (Decaying neutrino models have the added feature of a 
bump in the power at subgalactic scales to ensure early galaxy formation, a 
consequence of the large effective of the neutrinos before they decayed.) 

Although Fig. ^ gives a visual impression of which models are preferred, 
this can be put on a more quantitative basis. Fig. ^ shows trg is a sensitive 
function of n^. for CDM models with = 1, it is far too high at 1.2 for 
ns—1, but too low by ris ~ 0.76 with the "standard" gravity wave contribution 
{vt = Vs) or by rig ~ 0.60 if there is no tensor mode contribution. However, the 
shape constraint wants lower rig. In | [T^ , we marginalize likelihood functions 
determined with the COBE data (and smaller angle data) using a prior proba- 
bility requiring that F + z/s/2 be 0.22±0.08 and crs^^^r^ be 0.65tg J^ in order to 
condense the tendencies evident in Fig. ^ into single numbers with error bars. 
Threading the "eye of the needle" this way is so exacting that the error bars 
are too small to take too seriously. Sample numbers using only the 4-year dmr 
data and these priors for the 13 Gyr case are: = 0.76l g3 for h = 0.5 with 
gravity waves, 0.61^'q4 without, i.e., much tilt; for h = 0.7 and Vl^ac = 0.66, 
we get 0.99j;;g2 with and 0.94t Q4 without, i.e., little tilt; and when Hubble 
parameters in the range from 0.5 to 1 are marginalized over, the preferred in- 
dex is Us = 0.99l;o4 with gravity waves, 0.95j:;"o without, i.e., again little tilt. 
Combining sp9.l and .sk93+94 with dmr data does not change the values nor 
the error bars by much (it flattens a bit). When all smaller angle data is used 
plus LSS and dmr4 data, and Hq is marginalized, = 0.96 ±0.03 is obtained, 
again not very different from the LSS -I- dmr4 only result. In both cases, a 
nonzero ^vac is preferred. 

For the decaying neutrino model with = 1 to have erg > 0.5 we need 
F > 0.22, i.e., m^Td < 14 keV yr. The hot/cold hybrid model formula in 
eq. ( p^ ) is for one massive neutrino species. As Fig. ^ shows, an Hs—I h=0.5 
hot/cold hybrid model with fi^ < 0.3 would have erg > 0.8; however, even with 
a modest tilt to Ug — 0.95 this can drop to 0.7 for fli, — 0.25. (See also ref. 
.) That is, little tilt is required, in contrast to the CDM case. 
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Figure 6: This illustrates the accuracy and utility of the fitting formula for 
(Tg for a variety of inflation-inspired cosmological models. The lower figure also 
includes large scale structure and cluster abundance constraints. The symbols 
are defined in the text. 



It is also evident from Fig. g that the cluster data in combination with the 
dmr data stops h from becoming too high for a fixed age, but also would prefer 
a nonzero A value, with Hq 60 — 70 for 13 Gyr, and Hq 50 — 60 for 15 
Gyr. When the tilt is allowed to vary as well, the preferred values lower to 
very near 50 and 43, respectively, i.e., with little Qvac- h < 0.70 at 2a with 
gravity waves, h < 0.56 with no gravity waves for 13 Gyr; h < 0.56 at 2a with 
gravity waves for 15 Gyr. For the hot/cold models, the values near 50 and 43 
are preferred even more, even with very little tilt. 

For open CDM models, the COBE-determined as goes down with decreas- 
ing (and increasing h). These models are not so attractive because drops 
so precipitously with increasing h for fixed age {e.g., for 13 Gyr and Hq = 70, 
f^tot=0.055). Equation ([l^ ) has not been modified to treat open models (see 
e.g., (H). 

Texture and string models require a low as, < 0.5. At this stage, it is still 
unclear how much of a problem this is since non-Gaussian effects can modify 
the cluster distribution, partly compensating for the low as. 

Going to higher fc-bands, COBE-normalized spectra imply values for the 
"the redshift of galaxy formation" , of quasar formation, and the amount of gas 
in damped Lyman alpha systems etc. and this also significantly restricts the 
parameter space. However, it is more dependent on the role gasdynamics may 
play in defining the objects. 

If we are so bold as to assume that we now know the shape of the den- 
sity power spectrum over the large scale structure band, and the amplitude 
of the power spectrum on cluster-scales, then, in conjunction with the COBE 
anisotropy level, the range of inflation and dark matter models is considerably 
restricted. Whether the solution will be a simple variant on the CDM-|-inflation 
theme involving slight tilt (or more radical broken scale invariance), sta- 
ble ev-mass neutrinos, decaying (>kev) -neutrinos, vacuum energy, low Hq, high 
baryon fraction, mean curvature, or some combination, is still open, but can 
be decided as the observations tighten, and, in particular, as the noise in the Cf 
figure subsides, revealing the details of the Doppler peaks, a very happy future 
for those of us who wish to peer into the mechanism by which structure was 
generated in the Universe. 
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